Rabi oscillations in a quantum dot-cavity system coupled to a non-zero temperature 

phonon bath 
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We study a quantum dot strongly coupled to a single high-finesse optical microcavity mode. 
We use a rotating wave approximation method, commonly used in ion-laser interactions, tegether 
with the Lamb-Dicke approximation to obtain an analytic solution of this problem. The decay 
of Rabi oscillations because of the electron-phonon coupling are studied at arbitrary temperature 
and analytical expressions for the collapse and revival times are presented. Analyses without the 
rotating wave approximation are presented by means of investigating the energy spectrum. 

PACS numbers: 73.21. La, 42.50.Pq, 03.67.Lx 
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I. INTRODUCTION 



Semiconductor quantum dots (QD) have emerged as 
promising candidates for studying quantum optical phe- 
nomena 1]. In particular, cavity quantum electrodynam- 
ics (CQED) effects can be investigated using a single QD 
embedded inside a photonic nano-structure Q. One of 
the most fundamental systems in CQED is an atom in- 
teracting with a quantized field [3j, such a system has 
been an invaluable tool to understand quantum phenom- 
ena [H as well as for considerations on its applications to 
realize quantum information ||. Similar systems (in the 
sense of their treatment, applications, etc.) like trapped 
ions interacting with lasers @ have shown to be an al- 
ternative to develop techniques for quantum information 
processing 0] and the study of fundamental effects Q. 
Recent developments in semiconductor nano-technology 
have shown that excitons in QD constitute yet another 
two-level system for CQED considerations, and several 
successful experiments have been carried out @ , see also 
the referee [10(. Contrary to the atom-field interaction 
where dissipative effects can be fairly overlooked pro- 
vided the coupling strength between atom and field is 
sufficiently large compared to the dissipation rate due to 
cavity losses, the physics of a quantum dot microcavity is 
enriched by the presence of electro-electron and electron- 
phonon interactions. Thus, decoherence due to phonons 
may imply fundamental limitations to quantum infor- 
mation processing on quantum dot CQED [ll|. Here, 
we would like to analyze the effects of electron-phonon 
interactions on electron-hole-photon Rabi oscillations in 
cavity QED. 

As in Ref. 113. UM , we will not apply the Born-Markov 
approximation but will use a different technique 

for solving this problem. In particular we will make 
use of techniques commonly applied in ion-laser inter- 
actions. It relies on the rotating wave approximation 
RWA, and within this and the Lamb-Dicke approxima- 
tion the Hamiltonian becomes diagonal with respect to 
the phonon subsystem, which is shown in sec. |TTJ In the 
validity regime of the RWA, the decoherence effect on 



the inversion, due to the phonon bath, is analyzed in sec. 
IIIII and analytical expressions for the collapse and revival 
times are given. The zero temperature situation has been 
investigated in [13(, while here we study the effects due 
to zero as well as non-zero temperatures (causing the col- 
lapse of the revivals), and also how different phonon mode 
structures affect the decoherence. A different method to 
treat the non-zero temperature case was discussed in [l5| , 
where the collapse time is obtained numerically. The dy- 
namics beyond the rotating wave approximation, shortly 
studied here in sec. IIV1 becomes highly complex as can 
be seen from the energy spectrum. In what sense the 
phonon decoherence could be used as a possible resource 
for various applications is briefly mentioned in the con- 
cluding remarks. 



II. THE MODEL 

We assume a simple two-level model for the electronic 
degrees of freedom of the QD, consisting of its ground 
state \g) and the lowest energy electron-hole (exciton) 
state |e), with the Hamiltonian [l2} (%= 1) 



H = to eg a ee + cjara + g(a'&- + <r+6) 

+<3- ee ^ A fc (6j £ + b k ) + y^ j u k b\b k . 



(1) 



Here <r+ = |e)(<?|, a and b k are the annihilation opera- 
tors for the cavity mode and the fcth phonon mode, re- 
spectively. By transforming to a rotating frame, with 
frequency to 

Ti = Aa ee + g(a t <7_ + o+a) + & e ^\ k {P k + M+X/^jfrfc 

fe k 

(2) 

where A = oj eg — lu is the detuning. The transformation 

pea 



r = n< 



(3) 
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is used to obtain the Hamiltonian Ht — TliT^ 

H T = (A - A v ) a ee + ^2 Ukbffik (4) 

fe 

Ja-l[D b ( Vk ) + a + aY[Dl ( Vk ) 

k k 

with i]k = X k /u> k , A„ = X^fe^fe^fc IS the so-called polaron 
shift [HI and D~ b {"q k ) — exp (r] k & — Vjfij- For simplicity 
we look at the case A = A n to obtain 



where 



H T = g 



a^-IpS (Vk) + cT + aY[L>l ( Vk 



(5) 

Now by transforming to an interaction picture one ob- 
tains 



H, 



(6) 

The idea is now to apply the RWA on the above Hamil- 
tonian, equivalent to neglecting all rapidly oscillating 
terms. However, a closer look at (j6|) indicates that such 
a procedure is not straightforward. Normally uj k = ujQk, 
k = 1,2,3,... and expanding the above exponents give 
several time-independent cross terms leading to a com- 
plicated expression. However, it is know that the RWA 
is highly related to the Lamb-Dicke approximation [TtI ] . 
In the Lamb-Dicke approximation it is assumed that the 
variations of the exponents are smaller compared to the 
characteristic length of the phonon harmonic oscillators, 
typically rj k <C 1. One finds that in order to fulfill the 
RWA one normally needs to be in the Lamb-Dicke regime 
[13]. This observation helps us considerably in approxi- 
mating the above Hamiltonian, since now the cross time- 
dependent terms arising from different fc's vanish as we 
assume r\ k <C 1. In this case we can perform the RWA 
separately on each individual exponent to find 



Hrwa Q 



a' a. 



n La * fa*) + <5 + a n Lfi * 



(7) 



where L^ k are the Laguerre polynomials of order n k 
[lllii and g = ff exp(£/2) = <? exp (- £ fe 7^/2) is a 
rescaled Rabi vacuum frequency. We emphasize that to 
be consistent with the RWA and Lamb-Dicke approxima- 
tion, only terms up to rj^ should be considered when the 
Laguerre polynomials are expanded in the small param- 
eter rj 2 . This will be done in the following section. The 
parameter £ is sometimes referred to as the Huang- Rhys 
factor 20], and it is usually very small, ( « 1, [2ll |22| 
but it can become much larger, £ ~ 1, [23] and then 
our analysis would break down. The above equation is 
readily solvable, finding the evolution operator as 



U = U ee a e 



U eg (T- + U ge a + , 



(8) 



U ee (t;n k ) = cosfife )ft+ ii, 



U ge (t;n k ) = -tea 



n 



k.h 



U eg (t;n k ) = -iea 



I sin fl kt n + it 



Q k ,h 



+ 1 



and 



with 



Ugg(t;n k ) = COsOfe.ftt, 



(9) 
(10) 

(11) 

(12) 
(13) 



III. DYNAMICS 

Having the evolution operator, we can in principle cal- 
culate any properties we want, in particular we look at 
the Rabi oscillations for the two-levels system This quan- 
tity has as well been studied experimentally in quantum 
dot systems [22j]. By means of the inversion operator 
a z = \e){e\-\g){g\ 



W(t) = Tr{p(t)<r z } 



(14) 



where p(t) = Up(Q)W and p(0) is the initial density ma- 
trix, Rabi oscillations will be analyzed. This initial state 
is chosen as the fully separable state 



p(o)=n^(°)i°' e >< o < e i 



(15) 



with the QD excited, the cavity mode in vacuum and 
the phonon modes are all assumed to be in a thermal 
distribution 



pk{o)= KI 2 K)M = J2 (n +iw+i i nfc K nfe i- 

n fc =0 n fc =tr k ' 

(16) 



n k =0 



The inversion can be written explicitly as 



W ®= EIll^| 2cos [ 2 ~9H L n k ,(vl)t) , (17) 

{n k } k V fe' / 

where the summation goes from to oo and runs over 
all modes. For t = we have W(0) = 1 as expected, 
while for t ^ the sum will in general differ from 1. 
Interestingly we note that at zero temperature, T = 0, 
|c„J = 5 nk> o for all k and thus; at absolute zero temper- 
ature the system Rabi oscillations are intact but with a 
rescaled frequency [l3| . The expression for the inversion 
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can be simplified by taking into account that 77 <C 1 and 
we therefore expand the Laguerre polynomials 



W)=£ 



m— 



n \ if 



m J m \ 



1 — 7177 + 



n(n — 1) 



-n 



(18) 

Keeping only zeroth and first order terms in rj 2 , in the as- 
sumed Lamb-Dicke regime, the inversion can be written, 
after some algebra, as 



W(t) = cos(2gt)Ke JJ 



1 



rik + 1 — nue 



sm(2gt)lm ( J| 



1 



1 - n k c 



In the following we use dimensionless variables such that 
the quantum dot-cavity coupling g = 1, but keep it in the 
formulas for clarity. There are some special cases worth 
studying separately. 




FIG. 1: The atomic inversion W(t) as a function of time t 
for the single mode bath. The dimensionless parameters are 
■q = 0.01 and n = 40 (a) and 77 = 0.01 and n = 2 (b). 



TV identical phonon modes 



A. Single phonon mode 

In the simplest case consisting of a single phonon mode 
characterized by 77 and n, the inversion (|19|) simplifies to 



W(t) 



(n + 1) cos(2gt) - 77, cos [2gt(l - rj 2 )] 
(n + l) 2 + n 2 - 2n(77 -I- 1) cos(2gtrj 2 ) 



(20) 



The second term oscillates with a slightly shifted fre- 
quency compared to the first term causing a collapse of 
the inversion. When the two competing terms return 
back in phase, at times 



r W - k 



~ 9 ' 

m 



1,2,3, 



(21) 



inversion revivals occur. These arc perfect within the 
small ?7 2 expansion. For short times t, the inversion may 
be further approximated to give 



W(t) = 



cos(2gt) 



\ + %n{n + l)g" 2 r)H 2 ' 



Z~gtr? « 1, (22) 



and we conclude that the envelope function, determining 
the collapse time, is a Lorentzian with width 



1 



col 



(23) 



In fig. [T] we display two different examples of the atomic 
inversion (|17[) . The upper plot (a) has a large average 
number of phonons; n — 40 and 77 = 0.01. In the lower 
plot (B) the number of phonons is instead h = 2 and 
again 77 = 0.01. We can conclude that our results confirm 
that a lower temperature of the reservoir clearly increases 
the collapse time. 



By studying several identical phonon modes one may 
see the effect of multi modes in a simple analytic way. 
Instead of approaching eq. (fTT)) we go back to eq. ([5]) 
and let 77 = 77^ and lu = Lj k for all k. Let us introduce 
a N x N dimensional unitary operator connecting the 
boson operators b±, 62, 6 at with new ones Ci, £2, cjv 
such that 



1 N 



(24) 



k=l 



The unitary transformed Hamiltonian becomes 



H' T = g a)aJ) &1 (pVN) (r/Viv) +^c l t c i . 

(25) 

Thus, the problem relaxes to the single mode case with 
the scaled Lamb-Dicke parameter 77 — > r]\/~N. 



C. N different phonon modes 

In a more realistic model the phonon bath consists of 
non-identical modes, and depending on the model studied 
one has different spectral functions J(w) = J2k ^l^i^ ~ 
u>k) For the frequencies of interest here, J(w) has a 
simple power law behaviour [26[, resulting in a 
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J k > 



-2,-1,0,1,2,. 



(26) 



where k is a constant. The power s depends on mat- 
ter properties, for example; ohmic damping J(lo) oc lu, 
phonon damping J(u>) oc to 3 or impurity damping J(u>) oc 
lu 5 , but it also depends on system dimensions. We take 



4 



u>k = took, k = 1,2, 3, ... such that too determines the fre- 
quency spacing, while the thermal phonon distributions 
are determined from the average phonon numbers 





- 1 


exp (f ) 







(27) 



for some scaled temperature T. Often a frequency "cut- 
off" is introduced for the spectral function, but because 
the vacuum modes do not affect the system dynamics in 
our case, such a cut-off is not needed. Defining j k = 
8fi k (n k + l)g 2 i]l, and using the same arguments as above 
the collapse time in the small rj 2 expansion is given by 



n 



1 



lkt 2 co i 



)=2- 



(28) 



By further introducing 



r k = yj{n k + l) 2 + n\ - 2n k (n k + 1) cos(2gti]l), 
fc = arctan _ , , 



n k + 1 - n k cos(2gtrjl) 



(29) 



the approximated atomic inversion (|19p can be written 
as 



W(t) 



cos(2gt + 1 



(30) 



where 9 — J2 k and r = \\ k r k . In fig. [2j five examples 
of the atomic inversion are displayed for different powers 
s. Perfect revivals clearly occur for the s = 0, s = 1 and 
s = 2 cases, while for s = —2 and s = —1 the revivals 
are never fully perfect even at long times. This is be- 
cause if s = 0, 1, 2, ... we have rjj, — kuj^ — kwqF = rj^j 
where j is a positive integer, which does not hold if 
s = —2, —1. Note that the revivals are a consequence 
of the non Markovian treatment of the problem; the dy- 
namics is unitary and each phonon mode is considered at 
the same footing as the cavity mode and the QD them- 
selves. Therefore, also when a large number of different 
phonon modes are coupled to the QD-cavity system ex- 
hibits revivals, however less pronounced. 



IV. BEYOND THE ROTATING WAVE 
APPROXIMATION 

So far all the results have been derived within the RWA 
and the Lamb-Dicke approximation, which for grj/uj <C 1 
is expected to be justified. In this section we calculate 
the atomic inversion by numerically diagonalyze the trun- 
cated Hamiltonian of eq. ([5]) . The size of the Hamilto- 
nian is chosen such that convergence of its eigenstates is 
guaranteed. We restrict the analysis to the one phonon 
mode case, which already explains most effects. In the 
Fock state basis of the phonon mode, the matrix 



s=-2 



W(t) 




s=-2 



0.5 1 1.5 2 2.5 3x100000 




FIG. 2: Atomic inversion for different phonon modes obtained 
from the small rj 2 expansion (|19[) . The plots show respectively 
from top to bottom the cases with powers s = —2, —1, 0, 1, 2 
of eq. (I26|l . The other parameters are the same for all five 
cases; T = 0.2, k = 0.0004 and lu = 0.1. 



elements of the Hamiltonian are obtained by using the 
formula 



D mn = (m\D(rj)\ri) = 



e -N7V- m \P T L£- m (M 2 ), m < n, 



where L\ is an associated Laguerre polynomial. 



(31) 




0.5 



1 1.5 
(0 



2.5 



FIG. 3: The eigenvalues (|33[) as function of u, with r\ = 0.01. 
Dashed curves show fj,± and solid curves A±. 



For lowest order truncation of the Hamiltonian, one 
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keeps only a single Fock state of the phonon mode, and 
the RWA results given above are achieved [13] • In this ap- 
proximation, when the initial states of the phonon modes 
are on the form |c„J = 5 nkiVi i , the combined quantum 
dot-cavity system persists perfect Rabi oscillations with 
a rescaled vacuum Rabi frequency g — ► g Y[ n > L n ' (77^ ) . 
For zero temperature, all the modes are in the vacuum 
and the above approximation simply gives a rescaled fre- 
quency g — > g, which is the result presented in [13|. Thus, 
this simple derivation regains the results of [13| , and it is 
also clear how the approximation comes about and may be 
easily extended to non-zero temperature phonon baths. 
A deeper insight of the approximation (RWA) is gained 
by increasing the size of the truncated Hamiltonian. In 
other words, to go beyond the RWA one needs to in- 
clude more coupling terms arising due to the phonon 
bath. Considering an initial vacuum phonon mode and 
including first order corrections to the RWA, the Hamil- 
tonian can be written in matrix form (after an overall 
shift in energy) as 



LO 

~2 





Ao 


Ai " 





Ul 

2 


An 


D u 


Ao 


An 


U) 

~2 





An 


D u 





Ul 

2 - 



(32) 



with eigenvalues 

M ± = ^L±^ii ± y ( Z 3oo _D 11 -u,) 2 + 4Dg 1 . 

(33) 

These eigenvalues are shown in fig. [3] as function of ui 
for fixed r\; A± (solid lines) and /j± (dashed lines). The 
bare curves, e± = —Du±uj/2 and e± = Du±oj/2, where 
i = 0,1, are coupled by the non-RWA term Ai- This 
makes the crossings at u> = to be avoided, the degen- 
eracies are lifted. Thus, in the RWA, the difference in 
the values Ai for i — 0,1,2,... of the bare energies at 
u> = causes the collapse-revival structure. However, 
as the non-RWA terms start to dominate the spectrum, 
the bare energies are no longer proportional to ±u) and 
the full system dynamics will no longer show the nice 
collapse-revival structure. For which ui that this non lin- 
ear behaviour occurs, depends on g and r\. A rough es- 
timate can be derived by assuming that the linear de- 
pendence oiu/2 should overrule over the avoided cross- 
ing energy difference 2Ai> which after expanding in 77 
gives ui 3> Agrj. In other words, using our typical param- 
eter values of the previous section, we note that when 
u> < 1 the RWA is likely to break down. For small ui, the 
eigen energies are very densely spaced, and a small non- 
RWA correction will couple the different energies in an 



involved way. The number of states taken into account to 
correctly describe the dynamics is then growing rapidly, 
which can be seen in fig. 0] showing the numerically ob- 
tained eigenvalues for a 12 x 12 dimensional (a) and a 
22 x 22 dimensional (b) Hamiltonian. Clearly the more 
states taken into account, the more complicated energy 
spectrum. 




0.05 0.1 0.15 0.2 0.25 0.3 




0.05 0.1 0.15 0.2 0.25 0.3 
CO 



FIG. 4: The numerically obtained eigenvalues as function of 
lo, with rf — 0.05. In (a) the truncated Hamiltonian has size 
12 x 12 and in (b) 22 x 22. Due to the non-RWA coupling the 
curve crossings become avoided, which is, however, not seen 
in the plots of this sice. 

Finally, in fig. [5] we give one example of the numeri- 
cally obtained inversion in a regime where the RWA is not 
justified. Interestingly, for w <C 1 the state \g) (opposite 
of the initial state) of the QD is more strongly popu- 
lated, while for increasing uj the collapse-revival pattern 
appears and W{t) — > 1/2 in the collapse regions. 

V. CONCLUSIONS 

We have studied a quantum dot strongly coupled to a 
single high-finesse optical cavity mode by applying meth- 
ods usually applied in ion-laser interactions, namely the 
decomposition of the Glauber displacement operator in 
Laguerre polynomials. This allowed us to obtain results 
for the system when it is coupled to a phonon reservoir 
beyond the Born-Markov approximation. We have stud- 
ied several cases, including N identical and N different 
phonon modes, i.e. the case of non-zero temperature. Ex- 
pressions for the collapse and revival times of the Rabi 
oscillations have been derived analytically, valid in a large 
range of parameters. 

The analysis is carried out in the resonance condition 
A = A, ; , which in the RWA gives that no population is 
transferred between the QD and the phonon reservoir. 
If this resonance is not fulfilled, but other specific types 
of conditions are (blue or red detuned) , one may derive 
different types of effective Hamiltonians in the RWA, in 
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20 40 | 60 80 100 



FIG. 5: The inversion as function of w in the regime where 
the RWA breaks down. The revivals at t rev — n/rj 1 build up 
for increasing ui. Here rj 2 = 0.01 and n — 3. 



photon states will, just like in the case of the various 
phonon states, induce different Rabi frequencies affect- 
ing the collapse-revival pattern. Such situation may be 
of interest for state preparation or state measurement 
and is left for future considerations. 

We have also presented a short analysis of the dynam- 
ics in the non RWA regime. In this parameter range, 
the energy spectrum becomes very complex with cross- 
ing energy curves, and as the energy spacing between the 
curves is small for these parameters, many eigenstates of 
the Hamiltonian must be included in order to correctly 
describe the dynamics. Non the less, from such an ap- 
proach it is seen how the RWA results are obtained as 
a first order correction to the trivial situation, and how 
higher order terms cause avoided crossings between the 
energies, and they are therefore no longer proportional 
to ±a; around the crossings which will induce a "break- 
down" of the collapse-revival pattern. The higher order 
terms in such an expansion can be seen as a virtual ex- 
change of phonons [3(| • 



much resemblance with ion-trap cavity systems [3j . Then 
the effective models will be described by typical gener- 
alized two-mode Jaynes-Cummings Hamiltonians, which 
can often be solved analytically [29]. Another assump- 
tion made in the derivation is that the cavity mode is 
initially in vacuum, while for a general initial state each 
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